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We present an investigation of nonlinear interactions between Gravitational Radiation and mod- 
ified Alfven modes in astrophysical dusty plasmas. Assuming that stationary charged dust grains 
- - - form neutralizing background in an electron-ion-dust plasma, we obtain the three wave coupling 

I coefficients, and calculate the growth rates for parametrically coupled gravitational radiation and 

■ modified Alfven-Rao modes. The threshold value of the gravitational wave amplitude associated 

' with convective stabilization is particularly small if the gravitational frequency is close to twice the 

^ ^ ' modified Alfven wave-frequency. The implication of our results to astrophysical dusty plasmas is 

l- discussed. 
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• PACS numbers: 04.30.Nk, 52.35.Bj, 95.30.Sf 
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^ ■ I. INTRODUCTION 

cn '• 

^ ' There exists several mechanisms for conversion betwe en g ravitational waves (GWs) and electromagnetic waves 

plasmas. One of the most basic processes occurs 
when GWs propagate across an external magnetic field, which gives rise to a linear coupling to the electromagnetic 
field i, leading to the excitation of magnetohydrodynamic (MHD) waves in a plasma ^, Q. In order to excite 
^) perturbations with frequencies different from that of the GW, naturally nonlinear couplings must be considered. There 
\^ \ exist numerous examples of such mechanisms in plasmas, giving rise to, e.g. three-wave couplings between GWs and 
' electromagnetic waves. Wave coupling mechanisms involving GWs are studied for several different reasons. In some 
"q I cases, the emphasis is on the basic theory other works, the focus is on GW detectors flolfllL fl^ . on 

cosmology |l3l Il4l ^3 , or on astrophysical applications such as binary mergers |l6l | , gamma ray bursts |17|. pulsars 
^ ' ^3 supernovas 

In the present paper, we will consider gravitational wave propagation in plasmas containing charged dust particles 
ILJ ' |20|. The latter are prominent components in many astrophysical systems, and may contribute significantly to the 
. I dynamical properties of such systems . It has also been claimed that supernovae can be significant sources 

of dust particles [l^, although this claim is debated '25'|. Previous work involving dusty plasma-gravitational wave 
?-j interactions have considered general relativistic versions of the dust MHD equations ,2^^]. However, in cases 
where the dynamics is not dominated by the charged dust particles, other approximations are more useful |27j . In 
order to describe the modified Alfven mode (MAM) (or the Alfven-Rao mode 1271 ) > that can propagate in a magnetized 
dusty plasma, we will apply the infinite mass approximation for immobile charged dust grains. Thus, the only force felt 
by the charged dust particles will be the gravitational force, which is an appropriate approximation for a broad range 
of Alfven wave frequencies. Using the standard mode coupling theory |29| . we then obtain the coupling coefficient 
describing the nonlinear interaction between two MAMs and one GW. Using these results, we consider the parametric 
excitation of the Rao mode in the vicinity of binary mergers. Provided that the gravitational wave frequency is close 
to twice the Rao cut-off frequency the threshold value for the GW amplitude for the parametric interaction 
is much reduced compared to the ideal MHD theory assuming that the limiting amplitude is determined by 
non-dissipative stabilization. 
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II. BASIC EQUATIONS 



The metric describing a linearized gravitational wave of arbitrary polarization propagating in the z direction on a 
flat background can be written as 

ds^ = -dt^ + (1 + h+)dx'^ + (1 - h+)di/ + 2hy,dxdy + dz^ . (1) 
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A convenient choice is to introduce an orthonormal tetrad e^, where eg = dt and the spatial part is written as 
V — (ei, 62, 63) = Vq + Vg, where V and are given by 



Vg = -^{h+dx + hxdy,hy,dx - h+dy,0), 



(2) 
(3) 



and ei = (ei, 0, 0), 62 = (0, 62, 0), 63 = (0, 0, 63). 

Below we will derive the evolution equations for the matter and fields in the low- frequency MHD limit, in the 
presence of a gravitational wave. We consider a cold plasma composed of the electrons, ions and micron-sized charged 
dust particles. The momentum equation for the ions, when neglecting pressure and considering the electrons as 
massless fluid, reduces to 



{dt + V, • V) V, 



rriiC 



■(Vj - Ve) X B + gi 



(4) 



where 



is the gravitational acceleration of the ions 30]. For low phase speed (in comparison with the speed of light) the 
displacement current are neglected, and Maxwell's equations are written as 



a^B = VX (VeXB)-j^ 



(5) 
(6) 



where 



}E 

3b 



1 
2 

1 r 
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(E, - By)h+ + (Ey + Bx)hy 

{Ey + Bx)h+ — {Ex — By)h^ 



ei 



ei 



{Ey + Bx)h+ + {Ex - By)hy 



{Ex ~ By)h+ + {Ey + Bx)hy 



e2, 
e2. 



are effective currents due to the GWs, see e.g. Ref. 0. 

Combining Eqs. Q and (j^J and assuming the plasma to be quasineutral with negatively charged dust, that is 
qiUi — eue — qand = where qd > 0, we obtain 



{dt + Vi-\/)vi^ — - — (-ndqdVz + V x B - j^; ) x B -I- g,; 

en„mjC \ An / 



eneiriiC 

By using Eq. ((SJ we can eliminate Vg in © to obtain 



aB V X 



Finally, the ion continuity equation is 



V, - V X B H X B ^ - js. 

eUe AweUe eUe 



dtn^ + V • (n^Vj) = 0. 



(7) 



(8) 



(9) 



For later use it will be convenient to collect the gravitational terms on the RHS of the equations. Furthermore, we 
let the ions to be of charge qi — e. Equations iQ, © and ^ are then written as 

m.rieC ((9t V, • Vo) -f Udqd^i x B - -^(Vq x B) x B = -^(Vg x B) x B + riecmig, - x B, (10) 

An An 



dtB - Vo X 



n.iVi X B c(Vo X B) X B 



Aneue 



Vg X 



UiVi X B c(Vo X B) X B 



AneUe 



Vx 3b, (11) 



eUe 



and 



dtUi + Vo • {uiVi) = -Vg • (njVi), 



(12) 
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where the dust charge is defined by qd = Z^e. We note that in the absence of the GW-source terms, eqs. H10|) - 
(|12|l are MHD equations modified by the presence of infinitely heavy dust particles. Due to the dust charges in the 
quasi- neutrality condition, the relative contributions to the currents are modified as compared to ideal MHD. As a 
consequence a characteristic frequency fiij called the Rao cut-off frequency |0| enters in the linear wave modes. Wave 
modes propagating almost perpendicular to the external magnetic will be denoted as Alfven-Rao modes. In what 
follows we will be particularly interested in these modes, as they give the simplest way to intoduce a characteristic 
time scale (fi]^^), which is longer than the ion Larmor period, into the MHD equations. Furhtermore, as we will 
demonstrate, these modes are particularly easy to excite when the condition = 2riij is met, where ujg is the 
gravitational wave frequency. 



III. DERIVATIONS OF THE COUPLED MODE EQUATIONS 



In general a monochromatic wave of sufficient amplitude is subject to a number of instabilities, which can transfer 
the wave energy into other modes. If energy-momentum conservation allows for resonant three-wave interaction, 
typically this mechanism give raise to the most rapid instability |2fl| . Thus, in order to study this process, a system 
of three weakly interacting waves will be considered. Two dust MHD waves with frequencies and wavenumbers 
(wi_2,ki^2)j and a GW with arbitrary polarization propagating parallel to the background magnetic field with the 
frequency and the wavenumber (ajg,kg). We note here that a gravitational wave propagating in a finite angle to the 
magnetic field produces a linear coupling to the electromagnetic field 0. In the next step of the calculations, such 
linearly induced fields will complicate the description of the nonlinear mode coupling to a large extent. Thus our 
motive to let kg be parallel to the external magnetic field is to avoid a linear coupling between the GW and the MHD 
modes, and be able to focus on the nonlinear phenomena. As a prerequisite to obtain the nonlinear coupling, we 
first study the linear modes of the system (|l()|l - (|12|l . omitting the gravitational contributions. Considering the plasma 
waves to be plane wave perturbations on background quantities such that B = Bqz + B', n,; = no + f^' and = v', 
where So and no are constant and the primed quantities denote the perturbations, allows us to write the linear part 
of |Tni)-(ini as 



— I V X z — 



Ca [fc.B' - B^k] 



B' = 



r- [(k • v')z - fczv'] iniicCA kzk x B' 

ViTTnoniiCA h 



ey^Annomi uj 



— k • V , 



(13) 

(14) 
(15) 



where = ndqdBo/mi{no — Zdnd)c is the Rao cut-off frequency and Ca 
speed. The frequency matching is 

UJg = UJi +UJ2, 



7T,oBo/{no — Zdnd)\/4:nnomi is the Alfven 



(16) 



and since the gravitational dispersion relation reads ujg = ckg and Ca ^ c, the wavenumber matching can be 
approximated by 



ki 



ki 



(17) 



Thus we may consider the excitation of MHD wave-modes with wave- vectors that are almost perpendicular to the GW 
(or the external magnetic field). In particular we choose a coordinate system such that ky = and let |A:i^2z| ^ \ki.2x\, 
which allows the linear eigenmodes of the system to be represented by the following eigenvector 
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CAV4:-iTnomi 

The dispersion relation (DR) is now readily obtained and can be expressed as 

^2 ^ f^2^_^C'2^2^ 



(18) 



(19) 



which is the Alfven-Rao mode [23| , that reduces to the compressional Alfven (or fast magnetosonic) wave in the limit 



of zero dust-density. We note that the dispersion relation H19|l . together with Hlt)|l and H17|) . implies 



2u}i = 2uj2. 
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Next, assuming the dust MHD waves and the GWs to be plane waves with weakly varying amplitudes, we write 
/i+^x = /i+,x (O^'*'''''^"'^''*'' + c.c. and ip' — -0(i)e*''''^-^ ''~"^'^*-' + c.c, where c.c stands for complex conjugate, and 
tp' represents any component of B',v' and n' . Making use of the linear eigenvector H18|l as approximations in the 
nonlinear terms in the system (|10(I - (I12|1 . and keeping only the resonant part to second order in the amplitudes, we 
obtain the coupled mode equations 22j 

dtii,2x = «^ fl + ^) «2 ixK + 2nRvl^,hy, . (20) 



For the GWs we obtain, by using Einsteins equations linearized in h+,hx, and keeping only the resonant part of the 
energy-momentum tensor 



K 



4rj 



and 



dth+ = i—m^no ( 1 + — ^ ) vixV2x, (21) 



dthy = -4:K,m^nQ—^vixV2x, (22) 

UJ 

g 



for the + and x -polarization respectively. 

Noting that the gravitational wave energy density can be written as 



E, ^ W, + \hx |2) = ^ {\h+\' + \hx p) , (23) 

and the Alfven-Rao wave energy density as 

1 l-B P 

£^1,2 = Wi.2\vi,2x? = -^m^UQ {\,Vx.2xf + |wi,2yP) + = m,no | Wl,2:r P , (24) 

we can deduce three independent conservation laws from the coupled-mode equations For the total wave 

energy, i.e. 

dt 

the difference in Alfven-Rao wave quanta, i.e. 

d{Ni - N2) _ 



0, (25) 



dt 

and the sum of wave quanta, i.e. 

d{2Ng + iVi + iVa) = 
dt 



0, (26) 



(27) 



where we have introduced the number density of gravitational wave quanta Ng — Eg/Tiujg and of the Alfven-Rao wave 
quanta = Ei^2/fii^i,2- The existence of these conservation laws are equivalent to the Manley-Rowe relations [2^ . 

For simplicity, we have made the derivation of the coupled mode equations considering only time-dependent ampli- 
tudes. However, we note that a generalization to allow for weakly space-dependent amplitudes can be made by the 
simple substitution 9^ — > 9^ 4- Vg • V for each mode |2^, i.e. for a general slowly varying amplitude the coupled mode 
equations reads 

(dt + fgi,29,) ^1,2. = (^1 + ^ j + 2r!flw*_i,;ix , (28) 



{dt + cdz) h+ = i—m.ino { 1 + ) vixV2x, (29) 



and 



{dt + cdz) hx = -4:Kmino—^vixV2x, (30) 



where Vgi^2 — ki_2xC\/ lo i_2^ where ki^2x is the x-component of the wavevector for mode 1 and 2 respectively [3ll | 
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IV. SUMMARY AND DISCUSSION 



We have considered the nonlinear interaction between the modified Alfven (or the Alfven-Rao) mode and gravita- 
tional waves in a magnetized dusty plasma. In order to describe this process, dust MHD equations incorporating the 
effects of the gravitational waves have been derived. In particular, we have focused on the case where a gravitational 
wave of arbitrary polarization propagates parallel to the magnetic field. We have then calculated the three-wave 
coupling coefhcients for MHD waves propagating almost perpendicular to the magnetic field, in which case the latter 
wave modes obey the Alfven-Rao dispersion relation 1)19(1 . From the coupled mode equations we note that the GW can 
parametrically excite Alfven-Rao modes, which grow as exp{jt), where 7 depend on Qr, ujg, h+^x etc. Furthermore, 
it can be seen that 7 is roughly independent of the GW polarization, and of the order 7 ~ |/i+,x |wg. This estimate 
is the same as for the decay into high-frequency waves, see Ref or MHD waves, see e.g. Refs. 0, H^. The 
highest growth rate may be reached for approximately monochromatic gravitational waves, which could be produced 
by compact binaries close to merging, see e.g. 28], or during the black-hole ringdown [s^ ]. 

However, we note that the high gravitational amplitudes only exist during a limited time, and that the finite group 
velocity of the decay products has a stabilizing effect on the parametric process. To study this effect we consider the 
decay of a homogeneous intense gravitational wave which for definiteness has the polarization /i = . The amplitudes 
of the Alfven-Rao modes are assumed to have the form vi^ — vi^ e:xjp{iKx — iflt) and ~ i'2x^^Pi'''^^ ~ i^t)- 
Inserting this ansatz latter into Eq. ((28(1 . we immediately obtain the nonlinear dispersion relation 



{n-Vg,K){n-Vg2K) = --j:ll 



2 

9/1 I 



An 



16 V 



2 



2 



(31) 



Next, introducing Vg — \vgi \ — \vg2\, we deduce from that the growth rate of the Alfven-Rao modes is 




^-WK^vl (32) 

provided that the second term under the root sign is smaller than the first one. While Eq. ((32(1 formally shows that 
sufhciently long wavelength amplitude perturbations always are unstable, in reality there is a minimum wavenumber 
possible. This wavenumber i^niin is set by the shortest non-oscillatory scale of the problem, which may be either the 
inverse inhomogeneity scalelength 'Vno/no, i?~,rv or ictp)~^, where -Rcuiv is the background curvature and tp is the 
pulse duration time. Note that the pulse duration here refers to the time during which the pulse has a sufficiently 
constant frequency, such as not to break the frequency matching condition. The threshold value on the gravitational 
wave amplitude for the parametric excitation now becomes 

h+,,, = 8^^^_ (33) 

UJg 

Since the group velocity of the Alfven-Rao mode approaches zero when cu — > flij, we note that /i+trc — > when 
UJg 2ilii. Introducing the frequency mismatch 6uj = cOg — 2X1 and using Eq. 1(33(1 we find that the threshold value 
is 

/l+tre = -fi'minC'A. (34) 

Ug 

Thus, we conclude that provided the gravitational spectrum contains twice the Rao-frequency, parametric excitation of 
MHD waves occur more easily in a dusty plasma, that may exist in astrophysical systems, see e.g. Refs. |2lll23. 123ll2^ . 
as compared to a plasma without dust-particles. As a consequence of the low threshold, GW-induced Alfven-Rao 
modes may be excited at a comparatively large distance from the gravitational wave source. The aim of this study 
has been to shed further light on the gravitational-MHD interactions that take place in the vicinity of gravitational 
wave sources, such as collapsing compact binaries, quaking neutron stars, black holes during ringdown or supernovas 
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